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ABSTRACT 


The turning characteristics of the Bell 100 ton captured 
Air Bubble Surface Effect Ship (CAB/SES) are studied with 
the aid of the Oceanic Incorporated computer simulation. 
The craft's motions for straight runs, straight runs with 
waves and turning runs under calm water conditions, are 
compared with the final set of computer simulation turning 
maneuvers with waves involving changing speeds, wavelengths, 


wave amplitudes, and final rudder angles. 
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I. INTRODUCTION 


A. CAB/SES DEVELOPMENT 

The very latest design in displacement ships still 
manage to operate very inefficiently in the 40-knot-plus 
range. In fact the state-of-the-art in hull design peaked 
in the early part of this century. Ship designers have 
realized all along that ordinary displacement vessels 
sacrificed much efficiency for added thrust at speeds in 
the 40 knot range due to hydrodynamic drag. Propulsive 
thrust required increases approximately as the square of 
the desired velocity. To bridge this 40 knot gap, naval 
architects turned toward the surface effect ships (SES). 
Hovercraft, hydrofoil, air cushion vehicles (ACV), ground 
effects machines (GEM), and captured air bubble (CAB) craft 
all come under this generic grouping. A SES is basically 
any ocean going or amphibian craft which has the capability 
to operate in the air/water interface, thereby greatly 
reducing hydrodynamic drag. This study utilizes the CAB; 
in particular, the BELL 100B CAB/SES. 

The CAB was chosen over other SES due to projected 
state-of-the-art feasibility for high speed and sea keeping 
capabilities. Whereas hydrofoils and ACV's are size limited, 
the CAB is not. The CAB is mainly supported by an air 
bubble enclosed in the hard structure of a deck, two 
sidewalls, port and starboard, and flexible seals at the 


bow and stern. The air bubble is provided by a blower fan 
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system which provides enough pressure to support and lift 

the craft. Unlike the ACV cousin, the CAB retains its air 
pressure with the exception of leakage and occasional 

venting in higher sea states and maneuvers. The craft 

reacts much like a displacement craft when not supported on 
the bubble but maneuvering is extremely sluggish. Propulsion 
in several CAB/SES already built has been either waterjets 

or supercavitating propellers. 

When "on-the-bubble", there is a pocket or depression 
in the water equivalent to the gross displacement minus that 
part of the displacement which is being supported by the 
sidewalls and seals. As the craft increases speed and moves 
forward so too does the depression. This causes the same 
type of waves as a displacement hull. As acceleration 
continues, the CAB overrides the forward portion of the 
depression. This is called the secondary hump and a bow-up 
trim is obtained. Continuing acceleration causes the 
depression to slip aft under the stern seal so that a smooth 
flow results. Hump speed has been surpassed, wave drag and 
‘skin friction drag have been significantly diminished, and 
the drag versus velocity curve decreases to values well 
below drag experienced by displacement vessels. 

Over the past few years a joint research program has 
been carried out by the Navy and the Maritime Administration 
under the heading of Joint Surface Effects Ship Project 
Office (JSESPO). Presently that office is known as SESPO 


since the project is now sponsored by the Navy alone. 
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Several craft have since been designed, built, and tested. 
Among these are the XR-3, which is presently undergoing 
further research and evaluation at the Naval Postgraduate 
School in Monterey, California. Two 100 ton craft, the 
Aerojet General 100-A, propelled by a gas turbine water jet, 
and the Bell 100-B, a supercavitating propeller driven craft, 
are also undergoing evaluation by the SESPO. The ultimate 
projection of this program is to provide a CAB/SES in the 
multithousand ton class with cruise speeds from 80 to 100 
knots. 

The next step the project is expected to take is con- 
struction of a CAB/SES in the 2000 ton range. However, 
there are presently many questions such as seakeeping 
ability, safety, and dynamic behavior, which must be 
answered. Of the several methods for answering questions 
such as these, the most expensive would be to construct the 
full size prototype. On the other hand, towing tank model 
trials could be run, but probably the most effective method 
would be a computer simulation. In the long run the computer 
senuketiion data would have to be corroborated by empirical 
data. If the simulation were exceptionally good or at 
least a close approximation perhaps millions of dollars and 
many manhours could be saved. In the case of the CAB/SES 
the state-of-the-art rests in the two 100 ton craft already 
built. Instead of immediately building a 2000 ton craft 
based upon minimal data already obtained, Oceanics Incorpo- 


rated was contracted to provide a computer simulation, for 
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the loads and motion study, see computer program section for 
the Bell 100B,with the following overall objectives: 


1. To provide a computerized form of the six degree 
of freedom equations for SES craft that will 
yield time domain outputs of motions for the on- 
bubble mode of operation. 

2. To determine the motions of specified SES confi- 
gurations with special reference to human habit- 
ability, platform requirements (for military or 
commercial applications) and maneuvering require- 
ments for the on-bubble mode of operation. 

3. To establish a mathematical model for predicting 
maximum wave induced bending moments, both 
longitudinal and transverse, on SES craft to 
provide techniques that would yield slamming 
response for the on-bubble mode of operation. 

4, To determine the stability derivatives in the 
vertical plane from measured data of the 
vertical plane motions of existing SES craft 
in a random seaway. 


Direct comparison of simulated data with experimental data 
should then provide a reasonable measure of the accuracy 
of the simulation. The next logical step is to convert 


the 100B program to simulate a 2000 ton craft. 


B. BELL 100B CAB/SES PROGRAM 
1. Simulation Statistics 
‘All statistics given were held constant throughout 
Ree nan ber but may be modified. Parameters listed as 
variable were in fact varied throughout the procedure. 
Total Craft Weight: 209,999.3 lbs 
Center of Gravity: 


Longitudinal (ft fwd. Transom): 33.09 ft 
Vertical (ft above keel): 1«09 fT 


Kaplan, Paul, James Bentson, and Theodore P. Sargent, 

of Surface Effect Ship (SES Moti 

art I - Equations of Motion of SES Craft with Six Degrees of 
Freedom, p. ©, prepared for Department of the Navy surface — 
Effect Ships Project Office, Contract No. C-1-35055, August 1971. 
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Total Craft length: he eit 


Total Craft width: 33.8 ft (port and starboard 
sidewalls) 


Ventral Fins: 3 £t by 3 ft 

Plenum: 
Length waterline: 655315 Lt 
Length deck: 65. Sent 
Width waterline: SUG Dat 
Width deck: 268.00 £¢ 
Height average: Sieh 

Propulsion: Twin screws variable thrust 

Steering: Twin rudders 

ie Turn Rate: Variable 

Maximum: Variable 

Speeds: Variable 

Fans: 


Bowseal: 2 
Sternseal: 2 


2. Description of Program 

The Bell 100B computer simulation is a Fortran IV 
program primarily designed to calculate time histories of 
motions and loads of a surface effect ship, operating in 
the on-the-bubble mode in a seaway in six degrees of freedom. 
The program is of modular form in that each major element 
of the craft has its own subprogram, 18 in all. This feature 
enables the programmer to make a wide range of parametric 
changes by simply changing input data. Additionally, this 
modular format facilitates changes or additions to any 
part of the program. 

For any given program run, the adn program calls 
the input and initialization subroutine INCON. The program 


then calculates the motion time histories by calling the 
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various Meailar subroutines for as many time increments as 
specified to complete the test case. During the calcula- 
tions, the subroutines will, if requested, output variables 
selected by print option switches at desired print intervals. 

Reference (1), the users manual, contains a break- 
down on the contents of each subroutine. However, since 
the all important input data is handled by INCON, INCON will 
be discussed. INCON is composed of 19 separate blocks. 
Within each block there are one or more options or option 
control tags. Note in Table I, some of the options are 
listed as not required. The other options, those required, 
are considered to be the basic or specific data options 
without which the program will not run. 

The original program's output occurred on a print 
time basis dependent on which print option switches were 
utilized. Additionally in the original program, the user 
could request print options which produced vertical and 
lateral plane summaries. See the computer program below 
as modified for this study. 

3. Formulation of State Equations of Motion 

The equations of motion utilized by the simulation 
are formulated in terms of the various force or moment 
contributing elements. 

n 
F ne F, (1) 
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Table I. CAB/SES Incon Data Options 


BLOCK OPTION REQUIRED DESCRIPTION 


NUMBER 
001 01 Yes Timing information 
02 No Print options 
03 Yes Integrators and error tolerances 
04 No Summary prints 
05 No Motion option switches 
002 Ol Yes* Summary mass properties 
02 Yes* Discrete mass properties 
003 00 Yes Craft geometric input 
004 ol Yes Sidewall input 
02 Yes Appendage input 
005 00 Yes Sternseal input 
006 00 Yes Bowseal input 
O07 2> 3 OL. Yes Plenum geometry 
02 Yes Froude number input 
008 01 Yes Propulsion input 
02 No Engine out input 
009 ol Yes Rudder input 
02 No Rudder motion input 
010 00 Yes Aerodynamic input 
011 ol No Waves; amplitude vs. frequency 
02 No Waves; amplitude vs. wave length 
012 00 Yes Initial condition 
013 00 Yes End of case indicator 
014 00 Yes End of run indicator 
015 O1 No Shears and moments 
02 No Shears and moments 
03 No Shears and moments 
016 Not Used 
017 Not Used 
-018 00 No Case title card 
019 01 Yes Bowseal fans 
02 Yes Cushions fans 
03 Yes Sternseal fans 


%either but not both of these options 


16 


faa 1 
en 
gets Te 
ey i 
Yee 


motmqo ated) Yo dred Jon tad versie 


Fy ie "bubble, # "sidewalls, f Pow seal, + "rudder, 


+5 + F 
waves 


i 


propulsion, 83 "aerodynamic, 


+ F 


(2) 


i 


stern seal ei "gravitational, 


where F represents a moment or a force and i represents a 
particular component such as in the lateral direction or 
eee ered a 
The craft's trajectory is referred to an inertial 
reference frame, a right handed cartesian coordinate system. 


Thus, Xi, Vo? and Zo show the trajectory with the craft's 


fo) 
center of gravity at time t = 0 as the point (0,0,0). Roll 
(¢), pitch (6), and yaw (W) angles are also referenced to 
this inertial coordinate system. Displacements toward the 
bow are considered as in the positive Xo direction; displace- 
ments to starboard are considered as positive in the xo 
direction; aneclacements vertically down from the center of 
gravity are considered to be in the positive Zo direction; 
yaw (~) is positive with a turn to starboard; roll is 
positive with a heeling motion to starboard; and a bow up 
attitude is a positive pitch angle. Normal ship's motions 
including the CAB/SES have small values of pitch (@) and 

roll (¢). The main angle indicating the orientation of 

the craft relative to the inertial frame is yaw angle (Vp) 


which may have a large range if maneuveuring occurs. Motion 
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in the Zo aireceion, heave, will also be greatly limited. 
Since roll, pitch and heave values remain relatively limited 
in magnitude, overall complexity in simulating hydrostatic 
and hydrodynamic forces and moments is greatly reduced. 

The axis system describing the CAB/SES is assumed 
to translate with the craft through its trajectory. The 
craft's axis system is also a right handed cartesian coor- 
dinate system with its origin at the center of gravity. 
The vector relationship between the craft's axes and the 


inertial reference frame are given by: 


X, cos wp sin yp 0 X 

Yo = sin yp cos wp 0 YG 
Z 0 0 a Z 
fo) 


where the coordinate axis in capital letters indicate vector 
quantities in the respective directions. 

| Waves are also orientated with these two reference 
systems; the inertial frame to determine the wave properties 
and the craft's reference frame to determine the wave effect 
on the craft itself. Waves are assumed to be propagating 
in the direction defined by 8B, where 8 is the angle formed 
by the normal to the wave crests in the direction of propa- 
gation and the positive Xo axis. Waves propagate at an 


angle y relative to the ship where y is defined by: 


Y=B-¥ (4) 
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For example, if the craft were not maneuvering, Y would 
equal 8, since wy equals zero (0), and for waves coming 
from the bow y would equal 180 degrees. . 

Considering a linear velocity Vas a vector quantity 
with components in the x, y, and z directions, 


V, = iU + jv + kw (5) 
then, 
Py at gd BP Pe PL 
a 
= Fe m(iv + JV + kW) (6) 


Application of Newton's laws of motion with assumption of 
constant moments of inertia equal to rigid body values and 


small angle analysis leads to: 


Axial Force 


X= EF, = a (mU) - mVR + mWQ (7) 


r 


Lateral Force 


iG = 
Y= IF “Sar (mV) + mUR - mWP (8) 


Normal Force 


d_ 


at (mW) - mUQ + mVP (9) 


Z= rf. = 
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Rolling Moment 
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K = mM, = T¥xP (I, 
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+ (R° - Q°)1,) + (PR - Q)T, (10) 


Z y 


Pitching Moment 
y 


M = =M, Siae 7 Ty? CE - eee - (P + QR) I Ly 
2 2 : 
+ (P° -R )I ex + (PQ - R)ILe Ca) 
-- Yawing Moment... . : ee Seeaeain | 
: dh, e j . 
N = IN_ = ao sa rook + (Ty - T) PQ - (Q +t RP)IL 
2 2 = 
t(Q4 — P ) + (RQ - P)I ey (Pe) 
where 
h, Tx Taxy “Tee a 
h a) lod I = 
y yx yy yz : eS 
hy Sox ly lz R 


and U, V, and W are linear velocities along and P, Q, 
and R are angular velocities about the x, y, and z axes 
respectively. Also associated with these equations are the 


kinematic relations 


$ = P + wsine (14) 
ae cosé@ (15) 


¥ R+ @6sinht 
aaa aN 
¥ cos68 cos®@¢ (16) 
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Next, consider the nonlinear inertial coupling as indicated 
in Equations (10) and (11). When the above equations were 

in the process of being formulated for the computer simu- 
lation, a computational approach was carried out using 
analog computers to determine the significance of nonlinear 
inertial coupling using the nonlinear model of inertial 
coupling terms. Results showed negligible difference in 
responses and the decision was made to simplify the equations 
of on by Gore Ua ne ae oa eee Belew one ee WEE « 
Diigia for eoeammaleatar ame faaieeeea acne: program, ; 


Equations (7) through (12) may be simplified and written: 


U m 0 0 0 0 0 X | 
Vv 0 m 0 0 0 0 Y -— MUR 
W 0 0 m 0 0 0 Z 
= (7) 
P 0 0 0 Ty 0 ale K 
0 0 0 Oo. 0 M 
: yy 
R 0 0 0 oa 0 Le N 


The above form of the state equations allows establishing 
the equations of motion in state variable form as a set of 
first order differential equations. W, P, and Q are 


augmented by the simplified relationships from equations 


W=Z : (18) 
P=¢ (19) 
Q=6 (20) 
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thus introducing additional state variables. These states 
may then be transposed back to the inertial axes by use of 
Equations (3), and trajectories are determined through 


simple integration procedures. 
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II. PROBLEM STATEMENT 


The BELL 100B is assumed to be modelled by the computer 
program created by OCEANICS. The computer program is to be 
used to determine theoretically, the craft's turning 
characteristics in calm water and in simple waves. Turning 
characteristics relative to roll, pitch, yaw, and their 
rates of change are expected to compare favorably with 
available empirical data. 
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III. PROCEDURE AND DATA 


A. PHASE I DATA BASE 

Prior to development of the above problem statement and 
acquiring the necessary data, an initialized data base had 
to be established and trends, where trends existed, had to 
be noted. The data obtained in this study were obtained 
using the IBM 360/67 computer at the W. R. Church Computer 


Facility at the Naval Postgraduate School. 


By te Straight Runs | 
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bmimouEn tiils i of the initialization chase 4 
longitudinal speeds (U) of 40 knots through 105 knots siecle 
sive at 5 knot increments were used. Straight runs of 50 
seconds duration were made at each speed. Velocity was set 
at a constant value and thrust was allowed to vary as a 
function of surge only. For all speeds, except 95, 100, 
and105 knots, steady state values of thrust (T), pitch (86), 
and draft (d) were recorded and used as initialization values 
for a steady state run of 10 seconds. For the ten second 
run, all values remained constant; see Table II. 

Speeds of 95, 100, and 105 knots in actuality did 
not reach a steady state at the 50 second mark as did the 
other speeds. However, values of pitch, draft, and thrust 
were averaged over the last 20 second period. These values 
were in turn used in the 10 second run which again did not 
reach steady state. Thus, the values in Table II for these 


higher speeds are actually average values. Plenum pressure 
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Table II. Initialization Values 


ORDERED SPEED (KTS) ONE ENGINE THRUST (LB) DRAFT (FT) PITCH (DEG) 


40 5091.5 L202 23 
45 5490.2 1.223 .- 26 
50 6002.4 1.166 28 
55 6589.3 ~ 1.103 31 
60 7221.3 1.003 33 
°65 7872.8 ee -956 «35 
70 8521.0 874 -36 
Dace Ateoattgns. p2U IE Costern xqudte,ie + 058. 
80 | 9731.4 694 | 39 
85 10254.3 598 39 
90 10697.8 -498 -40 
95* 10839.1 - 362 ood 
100% 11646.8 337 -60 
105% 12631.5 ~322 -61 


*Average Values 
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was held at a constant 92.8 pounds per square inch for all 
runs. All succeeding simulation runs utilized these 
initialized values, thus forming a steady state reference 
point for system perturbations. 
2. Straight Runs with Simple Waves 

To further the data basis for turning maneuvers, it 
was decided to make straight trial runs with varying wave- 
‘length (A), wave amplitude (a), and speed (U). Initial runs 
were conducted with wave fronts coming from three different 
relative directions, B = 000° (astern waves), 8 = 090° (port 
beam waves), and B = 180° (bow waves). For Bach of these 
directions, speed (U) was varied (U = 40,50,60,70,85,90 and 
100 knots). Wave amplitude was set to a = 0.5 ft and wave- 
length A = 720 ft. Since this wavelength was ten craft 
lengths and the amplitude was relatively small, the situation 
approached that of an infinite wavelength. The initial 
conditions were set in accordance with data noted above. 

Table III shows the wave encounter period as 
experienced by the craft. The program computes initial 


encounter frequency in INCON as follows: 


= ae 
ree 2n(gr/2T) (u cosB)/A . (21) 
Negative encounter periods occur when the craft overtakes 
the waves in a following sea. For waves from the bow and 
the stern, no motion was noted in the lateral plane. Pitch 


angle (6) did oscillate in both cases at the approximate 
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Table III. Initial Wave Encounter Period 


SPEED (knots) BOW WAVES (sec) STERN WAVES (sec) BEAM WAVES (sec) 


40 5.613 -105.26 11.859 
50 4.960 - 30.342 11.859 
60 4.443 =a 20 11.859 
70 4.024 - 12.520 11.859 
85 3.525 = 3.691 11859 
90 3.385 - 8.887 | 11.859 
100 = _-:32136 - 6.656 11.859 
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wave encounter period with slight phase difference. These 
trends deviated slightly at U = 100 knots. At 100 knots, 
still no motion was noted in the lateral plane. With the 
encounter period reduced to 2.136 seconds, venting, presence 
of air entry gaps along the sidewall, was noted at the 
center of gravity by the fact that draft was shown to be 
negative at various, though seemingly non-frequency related, 
intervals. The sidewalls of the craft did not completely 
leave the water as noted by the fact that the plenum gauge 
peer did not fall below zero as was seen in later runs. 

Also noted in these runs where perderbatten forces 
acted only in the vertical plane, was the fact that the 
relative difference between the accelerations of bow and the 
center of gravity never varied by more than 0.05 G in the 
"worst case' runs of U = 100 knots. In spite of this small 
variance, the craft's attitude in the vertical plane, 6, 
rose and fell with the rise and fall of the wave encounter 
period with both bow and stern waves. 

For Wave ianent the port beam, the wave encounter 
frequency was the Same encounter frequency as would be 
experienced if the craft were stationary. The bow and stern 
waves caused the initial velocity of U = 100 knots, 168.89 
ft/sec, to decrease to U = 161.69 ft/sec and U = 166.41 
ft/sec respectively by time t = 25 seconds. However, the 
beam waves produced only a slight oscillation about a 
nominal drop in speed of about .4 ft/sec. Outputs for the 


beam wave case which vary significantly from the bow and 
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stern wave cases were pitch, roll, yaw rate, and displacement 
along the y-axis. It was interesting to note that an oscil- 
lation in pitch angle was produced when there were no 
explicit perturbation forces either in the z-direction or 
around the y-axis. The oscillation period appeared to be 
approximately that of the wave encounter period but slightly 
out of phase by one to two seconds. These oscillations were 


smaller than the oscillations noted in cases with bow and 


stern waves. 


“as _ Further Connlettion of data from the beam wave case. 
showed that draft measurements oscillated bcp pio snanediy 
the same frequency as wave encounter frequency but were 
180 degrees out of phase when considering maximum amplitudes.-: 
Roll angle (¢) and yaw rate (R) also showed a marked perio- 
dicity with respect to the encounter frequency. However, as 
speed (U) was increased this relationship to the encounter 
frequency seemed to diminish. Roll angle showed no discern- 
ible frequency whatsoever at 100 knots. 
The next ae taken was to reduce the wavelength of 
the simple wave. In the same manner as explained above, 
with seas from three directions, computer runs were made 
using a wavelength i’ = 360 feet or five craft lengths, with 
speeds U = 40, 50, 60, 70, 85, 90, and 100 knots. Results 
with the following exceptions were approximately the same 
as the above case: even at the 100 knot level the pitch 
angle showed a definite relationship to the wave encounter 


frequency oscillations; venting did occur for the bow waves 
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and the stern waves but the excursion in plenum gauge 
pressure was much greater in both cases; differences in 
bow and center of mass accelerations were greater; longi- 
tudinal velocity dropped off more in both the bow and stern 
seas case, but longitudinal velocity for the port beam seas 
show a slight but almost constant decrease, vice an oscil- 
lation over the 25 second problem period; maximum pitch 
excursions for the beam seas were about the same in both 
cases, +0.5 degrees, but this Seuneor was Wamanean 
less for the bow and stern sea. cases. | : , | 
The next step in setting the basis was to Rees c tne 
parameter variances over a wide range of representative 
values. To accomplish this, it was decided to utilize 
wavelength values of relative craft lengths. Thus, wave- 
length (X) was varied as one-half a craft length, A = 36 
feet, one craft length, A = 72 feet, two craft lengths, 
A} = 144 feet, and three craft lengths, X = 216 feet. In 
that the plenum average height is slightly over six feet 
and that drafts of 1.5 feet have been measured with only 
a 0.5 foot wave amplitude, wave amplitudes of a = 0.5, 1.0, 
and 2.0 fect were chosen as representative values. Further- 
more, it was noted that in the program, the input wave 
amplitude "a" was in fact an amplitude, but when considering 
effects of sea waves, the peak to peak value, the crest to 
trough value or the wave height is the measurement usually 
considered. In other words, the above values of wave 


amplitudes chosen have heights of 1, 2, and 4 feet 


30 


 aresite baw Woe: att 900 nb enter 116 beqann 
nate med tog ed? wat settsotay nt but tyne 
“alyeno ma sety .saaenseb dnatanon, teomia dud & 
dectq quaizen sholvsq asldetq baoase 2S ada 
‘Sigg mE seem ett tyeda wrew anda mand odave 
Ulireotitigts exw fotstuoxe atid dud vent 
25880 90 rede bas 2 
cl iilieseaira es 
avidainesstqe? Io same: eblw a a9v0 ® 


atitits o¢ bebloeb aw 72 ,aldd deldlamoses 


~ovaw ,nasd? = .cadtatiatl 2tane ewidteles. co: esoley 
of = % .ttyral Meee .s Dlarond es Bofray eaw 173 
;adtanet Nets ont gtewt Sf © ds .tdgaed Sist 

! nl fao* of5 © « ,eddpgewl Wess satis bas (25 
; tost «fe sove ¢ltdgtie ef ddgied egateva atinela, t 
| tine A3Siw cevvesem nped ovad geo? 2.1 20. estas gg) P 
.O.f .@.0 = & 30 sobpdiique evew  ebit tions *vaN ooh! 
~“wotinit .etvisw sv headmeterqet ae nénod6 etew sot 0.8 
yew Jat of? ,aageng ef? at dard beton saw ist. 


gifneilaton metw sad ,.@heifqus ns tos} al esw “s" obud 


hah 
a7 some sch) Goole Net ot Neeg odt-,aevaw ase Yo. edog] 


ome i oe WE a8 ” oy «2 Ota led evaw ed? to euvtey ducts 


vvow ‘) Meulov syode wit ,shsow tedvo ale 


5 
di 
- i 
= 
1” 


.& .L 36 eit hed ove nedodo 


| be ae 


respectively. The latter height leaves little room for 
heave mctions and therefore, when water contacted the top 
of the plenum, whenever fan horsepower went negative, or 
whenever the plenum gauge pressure went negative, the run 
was considered to be aborted. 

The same procedure, as followed for the above cases, 
where X = 720 feet and A = 360 feet, was utilized for the 
“varying wavelength and wave amplitudes. For each value of 
_ wave amp IT PUGe a = (0. Bigs das 0, ands?’ ae a total of 84 simu- 
lations. were run, varying wavelength, wave direction, and 
speed. An example of the data collected is shown for beam 
waves, in Table IV. Data for bow and stern waves were also 
collected, but the data were omitted due to later usage of 
beam waves as an initial condition for turns. 

Behavior of the CAB/SES in the above wave conditions 
was best compared by categorizing the waves by wave direction 
and not by wave amplitude. The most "well behaved" of all 
the wave directions was the beam wave. Well behaved in this 
case means that of all the wave directions, beam waves 
provided fewer run abortions than bow and stern waves. Of 
the 36 runs, documented in Table IV, only three runs were 
aborted due to water contact with the plenum deck. The bow 
and stern waves had twelve aborted runs each. 

For bow waves pitch angle maximum excursion was 
greatest for the longer wavelength, lower speeds, and higher 
wave amplitudes. The maximum pitch for bow waves was 


recorded at t = 20.5 seconds of the run using the wavelength 
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A = 216 feet and amplitude a = 2 feet at speed U = 40 knots. 
The frequency of oscillation of the pitch angle seemed to 
correspond roughly to the wave encounter frequency. However, 
in that print output was a one-half second interval, it was 
difficult ot determine phase delay times, and for shorter 
encounter periods it was impossible to determine any phase 
or frequency relationship at all. Starting with the one 
foot amplitude waves. and continuing through, no pattern 
could be determined with gauge pressure fluctuation for the 
bow waves men the ghee PLS SA CEP elon. that as speeds 
Pecreasea. ec aeed eane Rime Eie tone increased. This 
was also the case for stern waves. Maximum acceleration, 
pitch producing, differentials appeared between the center 
of gravity and the bow increasing with wavelength until 
4 = 144 feet and then dropped off considerably for smaller 
wavelengths. Interestingly final speeds of runs which were 
not aborted hit a low point at A = 144 feet. 

Stern waves, aside from producing the same number 
of aborted runs as bow waves, produced decidedly more 
adverse reactions in the craft. For all wave amplitude 
groupings, » = 144 seemed to affect the craft more adversely 
than any other wavelength. Maximum pitch and pitch excursion 
were produced with wave amplitudes of a = 0.5 and 1.0 feet 
and speeds of 40 knots. Additionally, the most violent 
plenum pressure gauge fluctuations, bow and center of mass 


accelerations differences, and speed loss were also noted 


35 


RaW 9E pLovredat brwose Hiadeone & biny pple 


nofsiuexs dotic bie dadiqowemisel d2unelevew tento” yaa 


0, i 


PS apt 


‘qodtode ‘tot brs ,nemtt ysfeb sasiig soncrege 9 
ennilg Vie ontateteb oF so rain haar ew we ahs 
ano 943 dviw gotsias2 Iisa ‘qidenoktaien | 
wisiteg on ,aquotds SasuH Lanes: bas bar dahoael 
ad? tot molzaytoull stiseesg egvay thw faieeclal r 
ebooq? es i214 aarsqooxs sidtzaod ort one a 


od whe sue yi 


sidT ,bowsetont aenolsautoult mum tntmamuatteam 


Nolienwleoss mumize Sevan ntsge tot onad “ee 
"edpen ef? owowsed barasqde aLeleneratith sarttsol nor 
fidgm dtansfeven ctjiw gntaeetant wod ef) Soe wank 
soflama 207% ¢«f[dascbleneo- Tis Beqgoth mets bin toot ¥ 
wiew doltdw env to abeogue Tank? vfgntsaezednt - ark 
,teat Mtg @ x om gatog wot a 3d” be. 

19Gann semua odd pclowheng movt eblas ,.sevaw 908 
sice yibehleet bsophortq ~eovaw wed aa arnt bad ‘£0 
ebytifana evew Lie tol .eteao oft at enotioset ° 


Ylearevbs o1ram fists sAY ¢oeTte 08 bomeom OL =. 4, 


jee? 0.1 tra ¢.0 = eo Yo setrdiiaona ovaw cidiw Seouboug 
ineloly J@Qm ed? -yifepottibbA © cegond 0% to abesqu 
2 


team to sot Sone wad , emehtentgoul? spies oni naaee nunetg 


yOeTJOrn OF LS Spe oof fhm f) Sria . peunexe?® fb “as 


lt ee eee 


for A} = 144 feet. As one would expect, there were no yawing 
or roll moments from either the bow or stern waves. Aborted 
runs due to negative fan horsepower seemed to occur whenever 
gauge pressures went above 120 pounds per square foot. 
Results for the beam waves are presented in Table IV. 
Generally sean. maximum pitch and pitch excursions 
increased with speed amplitude and decreasing wavelength with 
only two exceptions for non-aborted runs during the a = 2 
foot simulations. The frequency of oscillation seemed to 
_match wave encounter frequency, in. every case with little, or 
in some cases, seemingly no phase dur remende: Mecaverer ten 
variances between the bow and the center of mass were 
considerably less than values noted for either bow or stern 
wave cases, as were plenum gauge pressure fluctuation 
extremes. The low point for final speed loss in most cases 
appeared at values of i = 72 feet. The major difference 
between bow and stern waves ana beam waves, for the non- 
' maneuvering case, were the lateral component wave forces 
which produce roll and yaw moments. For the beam wave 
simulations, roll in every case decreased with speed. Also, 
values of roll increased up to id = 72 feet and then dropped 
off at X} = 36 feet for cases where speed was held constant. 
Both yaw rate and roll frequencies of oscillation were 
related to encounter frequency with varying phase time lag. 
The yawing effect was noted to be considerably higher for 
higher amplitude waves and then dropped off within amplitude 


groupings when moving from A = 72 feet to A = 36 feet. 
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B. PHASE II TURNING DATA 
1. Constant Flat Turns 

To provide a wide basis for turning data speeds of 
U = 40, 50, 60, 70, 85, 90 and 100 knots were each simulated 
BErrOuUrho burs Of Constant 5, 7.5, 10, 12.5, 15, 17.5, and 
20 degree rudder. For each rudder used under each speed the 
problem was allowed to run on a straight course for two 
seconds prior to moving the rudder. Initial condition values 
‘as determined above were not varied. At the two second mark, . 
the rudder was moved five degrees per second until the final 
rudder was reached. Thus, the final rudder angles were 
Beached at times 6 = 3.0, 3.5, 4.0; 4:5,-5.0,;.5.5,5 and 6.0 
seconds respectively for the above rudder angles. Figures 
1 through 9 show roll, pitch, and yaw for representative 
speeds U = 40, 85, ana 100 knots. In all cases the ordinate 
units were measured in degrees and the abscissa's units 
denoted problem time up to 25 seconds. To reduce confusion 
on the plots, only 5, 15, and 20 degree rudder values were 
plotted. For these turns it was noted that values of pitch 
and roll tended to oscillate about an average value. As 
time increased, the average value of roll/heel tended to 
decrease. This may be seen in Figures 1, 4, and 7, where 
the negative values of roll indicate a heeling motion to 
port for 5, 15, and 20 degree turns at HO, 85, and 100 knots. 
In the analysis of the yaw rate plot, Figure 10, yaw rate 


also oscillated around an average negative value, indicating 
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a turn to port. Yaw rate appears to be constant when 
considering the yaw versus time curves in Figures 3, 6, and 
9. 

In all cases, there was a steady decrease in speed 
throughout the 25 second turn problem; Table V. The more 
rudder used, the more final speed was down from the initial 
value. Average final roll angle, or heel, increased with 
“increasing speed; Table VI. Yaw and yaw rate reacted in 
like manner but also increased with rudder angle as well; 
see naples VII and VIII. Pitch angle was slightly more. 
erratic. For 40 and 50 knot runs, final average pitch 
angle for five degree turns were slightly less than the 
initial condition steady state pitch angles. In both cases 
as rudder angle was increased, pitch angle decreased. For 
all other speeds, the five degree final average pitch angle 
was higher than the initial condition value for the corres- 
ponding speed. As larger rudder angles were used the final 
average pitch angles increased and then dropped off at 
rudder angles of 7.5 degrees for 60 knots and 85 knots, and 
dropped off at 10.0 degrees rudder for 70 knots and 100 
knots; Table IX. Turning radii varied considerably for 
five degree rudders on a speed to speed basis. However, | 
turning radii seemed not to be dependent upon speed for 15 
degree and 20 degree rudders; Table X. 

2. The Roll Phenomenum of the 40 Knot Turns 
While plotting various turn data for the forty knot 


turns versus problem time, it was noted that values for roll 
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RUDDER ANGLE IN DEGREES 


15 


20 


5.0 


Uos 


10.0 


12.5 


15.0 


17.5 


20.0 


Table V. 


40 50 
66.7 ° 83:0 
60.7 74.4 
56.4 68.7 
Table VI. 
peo 
352 593 
-469 - 766 
541 -872 
tT -897 
- 580 871 
-547 -824 
-496 fST 


Final Time Simulation Speed (feet/second) 


SPEED IN KNOTS 


60 70 

99-1. (114.7 
87.1, 2 98.5 
197.0 S9L6 


Average Final 
SPEED IN KNOTS 


._. GOos 


-86 


1.10 


1.24 


1.25 


1.20 
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85 90 100 
¥36569 L43219 15352 
EV2250 196.49 12385 


101.0 104.2 109.8 


Roll Angle (degrees) 


e020) ORR, 100. 


CUChGe sey aleg7 


1.5958 2.03) 2.30 
2.19) 92-33) 2-50 
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RUDDER ANGLE IN DEGREES 


15 


20 


5.0 


7.5 


10.0 


2.5 


15.0 


17.5 


20.0 


Table VII. Final Yaw Angle (degrees) 


40 
16.4 
53.7 


70.6 


40 


SPEED IN KNOTS 

50 60 70 80 90 100 
20.5 122363 .026.35 .928.9 1429'.9 533.0 
64.8 75.3. 84.6 97.8 101.0 109.0 


S560 099.0 44954 ,129.6 194.5 145.5 


Table VIII. Average Yaw Rate (degrees/second) 


SPEED IN KNOTS 


BBP. 80, 0870580 7190., 4100 
8 | EGO) “algal cn aaa | 1.3 
iLAs) 1.5 1.7 1.8 1.8 2.0 
1.7 Delt tee <3 2.6 2.6 2.8 
Dee 2.6 2.9 3.3 353 355 
2.7 3k 13.5 3.9 4.0 4.3 
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RUDDER ANGLE IN DEGREES 


12.5 


15.0 


Initial 


15 


20 


Table IX. 


1130 = =1120 


X. 


40 50 
220 -279 
-210 -276 
.197 - 267 
-178 -250 
~154 -218 
25 177 
_ 2095 -126 
323 | 28 
Table 
40 50 
5800 5850 
1590 = 1600 


SPEED IN KNOTS 


60 


- 340 


«sol 


-350 


-328 


-298 


~252 


206 
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SPEED IN KNOTS 


60 
6000 
1620 


1140 


4} 


70 
.407 
440 
445 
428 
. 387 
326 
“253 


- 36 


70 
6100 
1640 


1130 


85 

- 603 
~652 
-590 
-579 
-526 


454 


1363 
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Average Turn Radii 


85 


6750 
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1130 


90 

646 
.728 
- 706 
-626 
-560 


-489 


393 


-40 


(feet) 


90 
6720 
1650 


1140 


Average Pitch Angle (degrees) 
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were not compatible with expected results; see Figure l. 

The value of roll for the first positive cycle, roll to 
starboard, for five degrees was greater than the value of 
the twenty degree turn by about 0.8 degrees. Furthermore, 
throughout the rest of the simulation, the amplitude of 

the five degree rudder oscillation was substantially greater 
than the amplitude for the twenty degree rudder. | 

To gather data to explain this apparent phenomenon, 
_initial results were first double checked both with smaller _ 
tolerances for the: Runge-Kutta-Merson numerical. integration 
technique and a constant step size; see computer program 
section. Results showed initial calculations to be correct. 
Next, values of roll were retabulated for rudder angles of 
five, fifteen and twenty degrees; Table XI. It was noted 
that the values for the fifteen degree rudder were the same 
point by point as the five degree rudder until time t = 3.5 
seconds; the rudder had stopped moving for the five degree 
turn at time t = 3 seconds and time 3.5 seconds the rudder 
was still swinging for the fifteen degree rudder simulation. 
Likewise, the values of roll for the twenty degree turn 
were the same at time t = 5.5 seconds as the values for 
the fifteen degree turn. 

The program was then modified, see the computer 
program section, to provide an optional output which included 
summations of axial forces (X), lateral forces (Y), vertical 
forces (Z), roll moments (K), pitch moments (M), and yaw 


moments (N); see Figures 11 through 16. Further modification 
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Table XI. Forty Knot Turn Roll Data 


RUDDER ANGLE IN DEGREES 


Five Fifteen 
0.0 0.0 
-09 09 
52 52 
1.07 1.18 
1.07 1.61 
47 1.53 
26 1.05 
58 -48 
27 0.0 
42 13 
1.02 18 
1.15 73 
74 1.16 
14 1.18 
25 84 
20 42 
.20 23 
68 35 
91 64 
77 89 
39 91 
05 74 
05 sD2 
13 42 
45 48 
69 -63 
-70 74 
51 75 
25 -66 
10 -56 
10 aby, 
32 55 
2 62 
61 - 66 
54 66 
. 38 61 


43 


San 


SASAkAs 


- 


ea 


was made to the basic program to provide the components of 
the roll moment (K). These include the bowseal, the stern- 
seal, the rudder, the side walls and appendages, the pro- 
pellers, the bubble, as well as components from wave and 
aerodynamic inputs. As above, with roll angles, the summa- 
tion of moments about the x-axis were equal through time 
3.0 seconds. At t = 3.5 seconds, the summation of monene 
for the five degree turn case jumped to over 44,000 ft. 1b. 


where as the summation of moments was only about 14,000 ft. 


‘- lb. for the twenty degree turn. See Table XII for precise 


tabulation results of times t = 2 seconds through t = 7 
seconds for the summation of moments about the x-axis. 
3. Turns with Simple Waves 

The final steps of this procedure included deter- 
mining what and how many simulation runs to make to insure 
valid data. Original intentions were to establish a field 
of various wave lengths, speeds, and wave amplitudes from 
all three wave directions as utilized in the above procedure. 
However, during those runs it was noted that the craft's 
reactions to parametric changes were highly dependent upon 
wave encounter frequency. The craft reacted in a much 
more seaworthy manner with beam waves. For higher veloci- 
ties and amplitudes and smaller wavelengths, simulation 
runs were often aborted due to loss of plenum pressure, 
water contact with the deck of the plenum or due to negative 
plenum fan horsepower caused by high pressures within the 


plenum. It was noted that the craft will experience 
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Table XII. Summation of Principle Roll Producing Moments 


Five Degree Rudder Simulation 


M M M M RUDDER 
x rudder bubble sidewall 

(£t.1b.) “Gee stb.) (ft.1b.) (ft.1b.) (deg) 
0.0 0.0 0.0 0.0 0.0 
32,522.0 42,011.0 1700.0 10,359. 25. 
19,994.0 63,742. 4872. 48,782. 5.0 
44 ,603.0 36,589. 20132. 91,491. 5.0 
49,689.0 28,959. 20510 89,133. 5.0 
10,847.0 29,719. 9673. 45,510. 5.0 
29,839.0 SZ S491. 3517. 854. 5.0 
45,134.0 35), 27. 4501. 14,864. 5.0 
28,444 .0 37,465. 4676. 6,630. 5.0 
5,060.0 39,503. 7010. 48,214. 5.0 
32,804.0 40,868. WI D22 = 82,627. 5.0 


Fifteen Degree Rudder Simulation 


SAME AS ABOVE 0.0 
2.5 

5.0 

14,672. T8594 21,788. 104,762 USS) 
37,159: 95,386. 29,625. 147,688 10.0 
26,000. “215,132. 28,287. 158,146 12.5 
1,588. 138,015. 19,880. 146,565. 15.0 

1,414. 120,520. 10,047. 124,223. 15.0 

21,690. . 4 E19, 312. 1,469. 98,770 15.0 
30,961. 121,867. 1,458. 90,146. 15.0 
18,247. 126,382). 3,159. 107,645. 15.0 


_ Twenty Degree Rudder Simulation 


SAME AS ABOVE 0.0 
2.5 

5.0 

7.5 

10.0 

12.5 

1,580. 162,916. 15.0 
28 ,056 187,987. 11,594 140,705. 17.5 
31,414. E705 725i 9,062. 101,486 20.0 
25 I22 6 168,321. 13,316. 190,036. 20.0 
21,524. 170,424. 14,910. 197,390. 20.0 


NOTE: All moment producing elements are not included. 
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varying encounter frequencies with any initial wave direction 
utilized since wave encounter frequency is a function of 
wavelength, longitudinal velocity, and relative wave 
direction. For each simulation where the ship is maneuver- 
ing an instantaneous encounter frequency is determined. 
Thus, it was decided to start each problem with waves from 
the port beam, allow the craft to follow a straight course 
for the first two seconds and then put the rudder over at 
time t = 2 seconds with a rudder swing of five degrees per 
second as was done above for the flat turns. Again initial- 
7 ized data as determined in the ns part of this procedure 
were used. ¥ + | . | 
At this point another cecision had to be made in 
order to provide enough useable data. The original plan was 
to make sets of twelve computer runs for each of the speeds, 
U = 40, 50, 60, 70, 85, 90, and 100 knots. Six-hundred- 
seventy two simulations averaging seven minutes each were 
nearly out of the question. On a second cut, it was 
decided to use each of the above speeds with eight runs. 
Each run consisted of seven cases with varying rudder angles 
of five degrees through twenty degrees in 2.5 degrees 
increments. The eight runs used various combinations of 
previously used wave amplitudes a = 0.5, 1.0, and 2.0 feet 
and wavelengths \ = 216, 144, 72, and 36 feet. For each 
eroup of runs, A = 144, a= 2: } = 72, a= 23.A = 144, a= ls 
and } = 72, a = 1, was deleted. These values were deleted 


not only because of the computer time requirement, but also 
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due to the fact that these four sets provided repetitive 
values of wave amplitude to wavelength ratios: 


Max slope = 360°° = (22) 


as determined in the subroutine INCON. 

Table XIIIT is presented to show the format for data 
collected at all speeds. Particular turning maxima and 
Manima are noted in Table xiv. The terms Ad and AQ in 
Table XIV are maximum roll and pitch excursions respectively. 
-Maximum values were absolute values while the minimum values — 
were the minimum of the maximum values of collected data. 

In general, roll seemed to increase with decreasing wave- 
length until A = 36 feet. At this point roll values 

decreased significantly. In all cases, roll values were 

less for wavelengths of 36 feet than they were for wavelengths 
of 216 feet. 

Pitch generally decreased from i} = 216 feet with 
decreasing wave Lerieth®., Side slip angle, the angle formed 
by the direction of force summation and the ship's longitu- 
dinal axis, varied with increasing rudder angle. Rudder 
angles of five degrees, fifteen degrees, and twenty degrees 
corresponded to average side slip angles of 2.3, 5.8 and 
7.3 degrees respectively. Yaw angle and yaw rate increased 
with increasing rudder angle and speed. Further observations 
not directly associated with turning data were gauge pressure 


readings and final speed relationships. Again, as with 
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Table XIV. Maxima/Minima Summary for Turns with Waves 


MAXIMA MINIMA 
Value Simulation Value Simulation 
(deg) U a x Srud (deg) U a X Srud 
(ft/sec) (ft) (ft) (deg) (ft/sec) (ft) (ft) (deg) 
ROLL 4.84 40.0 0.5 72.0 5.0 1.68 40.0 0.5 36.0 5.0 
ROLL 9.35 40.0 0.5 72.0 5.0 0.90 100.0 0.5 36.0 5.0 
PITCH 3.09 40.0 1.0 216.0 5.0 0.66 50.0 0.5 36.0 5.0 


PITCH 4.51 50.0 1.0 216.0 15.0 0.27 100.0 0.5 36.0 15.0 


opt 


mS £ a UD (aod) tox A a Os 
(gob) (23) (33) (owals2) gob) C93) (32) Caen ae. 


0.2 0,88 2,0 0.08 a8.£ 0.2 0.89 2.0 008 8a 
et, hl OR O,0E 2.0 0.002 08.0 0.2 o.st 2.0 0a 20.8 
7 2 Oat ; 2.0 0.02 88.0 0.8 OMe wa 9.08 ROE b | 
O.8L 0.ar 2.6002 NS.0 0.8L 0,846 Ot 0.08 te, a nN 9 
a f ~~ 


- 


Le 


straight runs with waves, runs were aborted when plenum 
gauge pressure surpassed 120 psf. Runs with wavelengths of 
36 feet provided a slight decrease in gauge pressure. Speed 
drop-off was also noted with increasing time, speed, wave 
amplitudes and decreasing wavelengths with the exception of 


36 foot wavelengths where final speed values increased. 
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IV. DISCUSSION/EVALUATION OF DATA AND PROCEDURE 


The original plan of attack for obtaining turning data 
was both broad and sound considering the fact that no 
information other than Reference 1, the user's manual, was 
available for simulation information. The data base was 
established by determining relationships for straight runs 
with simple waves, and by determining craft reaction to | 
calm water condition turns. The procedural part of the 
original plan was to utilize a wide set of ae points by, 
varying speeds from 40 knots to 105 knots, varying rudder 
angles from five degrees to twenty degrees, varying wave- 
lengths from one half of a craft length to ten craft lengths, 
varying wave amplitudes from 0.5 to five feet, and varying 
wave directions from 180, 90, and 000 degrees relative. 

The sheer magnitude of possible data accumulation demanded 
that data points Be minimized and chosen so as to provide 

a wide spectrum of relevant data. Thus, after initial 
conditions had bean found for speeds in five knot increments 
from 40 knots to 105 knots, the field was cut to only seven 
speeds for the straight runs with waves. Prior to proceeding 
to the turns with waves, it was decided to limit initial wave 
direction to port beam waves. The beam wave data from the 
flat turns were further subdivided and only 40, 60, and 

100 knot information was extracted. Finally, of 84 possible 
Simulation runs for each of the seven speeds for the turns 


with waves simulation, 28 of these were deleted and only a 
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total of 24 runs, 5, 15, and 20 degree rudder turns, were 
recorded for the broad analysis. In the final analysis, 


only the 40 knot turns were presented in this study. 


A. PHASE I DATA BASE 
1. Straight Runs 
The straight simulation runs with no waves set the 
initial conditions for all succeeding runs. Although steady 
‘state conditions were not determined for the higher speeds, 
U = 95, 100, and 105 knots, it was decided to accept the 
average values as steady state conditions. 
‘2. Straight Runs with Waves 
The first simulations with waves utilized long 
wavelengths and small wave amplitudes with the intent of 
extracting information due to nominal perturbations of the 
wave effects at what approximated a wave of infinite length 
relative to the craft's dimensions. Unfortunately, the 
first wavelength used, A = 720 feet, yielded ‘more reaction 
than had been expected. Rather than Mmeresiime wavelength 
to approach the originally desired relative infinite wave- 
Hereth, it was decided to halve the wavelength to five craft 


lengths or iX = 360 feet. The final decision to use wave- 


lengths of A = 216, 144, 72, and 36 feet was based on two 
main ideas. Firstly, there had to be a choice of wavelengths 
that would provide relevant data. Secondly, actual seas 
encountered by such a craft would not usually be of the 


extreme wavelengths as encountered in open seas by larger 


vessels. 
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In the procedure section above, it was moved that 
wave encounter frequency played an important role in the 
ship's stability characteristics. Roll, pitch, and yaw 
all oscillated with some relationship to this wave encounter 
frequency. The relatively gross one-half second print inter- 
val used to determine the values for the beam wave simula- 
tions showed only a general relationship and negligible 
phase information. Figure 17 shows wave amplitude for a 
straight 40 knot run with A = 216 and a = 0.5 feet. When 


wave amplitude was compared with Figures NO tanid! NOt lots 


*..0f roll and pitch, it was. readily seen that. roll and pitch 


were definitely wave encounter frequency related with only 
a slight delay time. Interestingly, the roll and pitch 
producing moment plots, Figures 20 and 21, showed no such 
clarity for frequency or phase relationships due to the 
irregular formation of the moment curves. However, when: 
the roll and pitch curves were compared to their respective 
moment curves, it was noted that they were approximately 
180 degrees or about one-half of a wave encounter period 
out of phase. 

The loss of speed noted in the wave procedure was 
due to the increased drag caused by the wave. Figure 22 
shows wave drag in the x-direction with varying rudder 
applications under calm water conditions. The important 
part of Figure 22 for this example is that it illustrated 
an initial wave force which was constant until the rudder 


was put in motion at the two second time of simulation. 
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It is this initial wave drag which when coupled with the 
additional hydrostatic drag forces caused loss of speed. 
As velocity decreased more wave drag was experienced, see 
reference [2], and speed was further decreased. The end 
result was that eventually a steady state would be reached 
whereby the initial thrust as set into the simulation for 
the desired initial speed would eventually equal the total 
drag forces thereby reaching a steady state speed. Times 
for bhis steady state SE QUE NO were Genet ce nabry longer 
than the 25 second simulation runs used in this report. 

For the bow and stern waves the roughest perturba- 
tions experienced by the craft were simulations with A = 144 
feet. Beam waves similarly provided maximum distrubances 
at X41 = 72 feet. The cause of this had to be attributed to 
the craft's dimensions. For longer wave lengths the craft 
experienced a rise and fall with the wave encounter fre- 
quency. However, for wavelengths of equal size to or smaller 
than the dimension in question, the 75 foot length for the 
bow and stern waves and the 36 foot width for beam waves, 
angular displacement excursions were greatly reduced due 
to the rapidly changing wave forces. The fact that these 
angular excursions were less was not in itself a fortuitous 
occurrance. As noted in Table IV, the only aborted runs 
were at the lower, 36 foot, wavelengths for the beam waves. 
Additionally, all of the aborted runs for bow and stern 
waves, not shown, occurred at wavelengths of 72 feet or 


less. The reason for these aborted runs was the excessive 
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angular accelerations experienced by the craft due to 
increased wave encounter frequencies and increased wave 


amplitudes. 


B. PHASE II TURNING DATA 
1. Calm Condition Turns 

Most of the data collected did not indicate any 
unexpected results. Roll, pitch, and yaw rate all exper- 
ienced some type of decaying oscillation commencing with 
Same) t1=,.2 seconds, the beginning of rudder motion.” This 
_ rudder motion was considered Ee be i, ramp TPB, me he 
reached its final desired aul and remeiied a Secbeae 
at that value throughout the simulation. With the cessation 
of rudder motion, force values due to rudder change also 
began an oscillatory decay. Roll characteristics will be 
discussed in greater detail below, but for the most part 
roll attempted to reach a steady state heeling attitude to 
port. Pitch final values, as noted in Table IX and illus- 
trated in Figures 2, 5, and 8, for 5, 15, and 20 degree 
turns with speeds of 40, 85, 100 knots, oscillated about 
values less than the initial condition values. Figures 2, 
5, and 8 indicate that the more rudder applied the greater 
the average pitch angle. Yaw as opposed to roll and pitch 
seemed to be fairly steady when viewing Figures 3, 6, and 
9. However, Figure 10 shows that a 40 knot simulation did 
experience oscillations in yaw rate producing more oscillation 


with the smaller rudder angles. 
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2. The Roll Phenomenom of the 40 Knot Turns 
As noted in the procedure, positive peak overshoot, 

roll to starboard, was much greater for five degree rudder 
than for either 15 or 20 degree rudders; Figure 1. The 
whole problem was eventually resolved by calculating the 
breakdown of the roll producing or the components.of the 
moments about the x-axis. When data such as found in Tables 
“XL ana Tor pene properly correlated on a time basis, it 

was seen that values of the components were exactly equal 
Pet) immediately eee ene, euddes Sr enpedu its aa Thus’: 
‘at times greater. than three seconds for the five degree 
turn, and times greater than five seconds for the 15 degree 
turn, a different problem was generated than for the 15 
degree and 20 degree turns respectively. In that the clock- 
wise motion of the rudder produced greater negative roll 
moments with increasing rudder at print time 3.5 seconds, 
the negative moment component due to the rudder was much 
Breater for the 15 and 20 degree rudders than for the five 
degree rudder. Si iayiy, at print time 5.5 seconds, the 
negative moment for the 20 degree rudder was significantly 
larger than the moment for the 15 degree rudder. As the 
craft begins to turn the dynamics are such that the summa- 
tion of forces, and for that matter the direction of motion, 
is not always directed along the craft's longitudinal axis. 
In fact, for a turning maneuver with constant rudder a new 
line of force is generated. This direction of motion 


relative to the craft's heading is known as side slip angle. 
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In the simulation program, side slip angle is simply 


determined by 


aI< 


Side Slip Angle = C23) 


in radians. The initial motion of the rudder provides this 


; Baan in. SHES COSI eh Sie At eine point, from BAL ‘XIII, 


it was noted that with shamuaceenas sideslip. angle ie be ive: Cte 


_ component of roll moment due to the sidewalls increased 
until the moment component peaked within spceiaale of a 
second after the rudder motion had stopped:- Average side- 
slip angles were 2.1, 5.8, and 7.2 degrees for 5, 15, and 
20 degree rudder turns respectively. Obviously, the in- 
creasing side slip angle produced increased moments with 
higher rudders, however; as noted above, these values did 
not peak until after cessation of rudder motion. Thus, at 
time 3.5 seconds the slightly large positive moment due to 
the sidewalls coupled with the notably smaller negative 
moments for the rudder produced a much more positive roll 
reaction, roll to starboard, than in the cases of the 15 
and 20 degree turns. 
3. Turns with Simple Waves 

The results of the first runs of turns with waves 
showed that the best conditions generated by the beam wave 
case soon deteriorated as the craft's heading turned toward 
the wave front. Due to the one-half second print interval 


decided upon for all the runs, frequency and phase analysis 
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for turning characteristics were not carried out. Figure 
23 plots the wave amplitude for a 40 knot 15 degree rudder 
wunewith A = 216 ft and a = 0.5 ft. Notice in particular 
how encounter frequency had increased near the end of the 25 
second run. Figures 24 through 26 show the same 40 knot 
turning run compared to a 40 knot turn without waves in 

roll paved ehig in Soe rene nt: Notice that the frequencies 
ee the Geeny ine beelliations mencned™ cor ‘Some. , déeree “the: eee 
ee one ne Ey aueMctee Tor no waves with only Slight 
time delay. Notice etse that in Figure AS. the turning 
_ path for the wave-.case seems to have a much smaller radius 
than the radius for the no wave case. This is attributed 
to the z-axis moment added by the wave. Although cases 
were not run past the 25 second mark, it is assumed that 
the turning radius for the wave case would be correspondingly 
larger when the z-axis contributed by the wave acts against 
the direction of turn when completing a 360 degree turn. 

As with straight run waves, turning runs with waves 
seemed very Gepenacne on wavelength. Slower runs with long 
wavelengths yielded maximum perturbations in roll, pitch, 
and roll rate. For turns with waves, the drop off point 
for some parameter excursions, and the minimum for others 
was for wavelengths of A = 36 ft. This seems reasonable 
since as in the above straight run simulations with beam 
waves, the beam dimension of 36 ft produced the same results 


when considering wavelengths. 
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V. CONCLUSIONS 


A. TURNING CHARACTERISTICS 

The study of any craft's turning characteristics in 
the final analysis must be based upon actual turning data 
generated by a prototype craft. For the purposes of this 
~ gtudy, corrélating data could: have: generated more ‘definite °.:-°. 
conclusions as to the validity of turning characteristics 
presented above. "However, since no such data has been made 
available only generalizations can be made as to the 
simulation capability of the Oceanic's Seoeeeae 

Basically, turning characteristics seemed to be compatible 
with previous knowledge and past experience. From a controls 
point of view, data tends to indicate that with this 
particular CAB/SES as modelled would be optimally utilized 
at higher speeds using larger rudder angles. Of course, the 
dangers of high speed maneuvering are obvious when considering 
increased and varying wave encounter frequencies. In that 
the craft was deshoned for high speed operations, control 
efforts should be directed at reducing both the acceleration 
and the maximum excursions of attitude parameters. As it 
is, the simulation showed only such small changes in these 
parameters that if the same perturbations were experienced 
by a displacement vessel, they would in general go unnoticed. 
However, since the CAB/SES is not a displacement vessel, 


these stability parameters are much more important. 
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overatl the simulation seems to be accurate when 
considering horizontal plane trajectories. As far as the 
above simulation runs were concerned, data obtained for 
vertical plane motion, heave in particular, leave room for 
some doubt. Original runs used the variable Runge-Kutta- 
Mersen numerical integration. As aves were introduced 
to the Samulatior ens) more eee vonae was Seta re due 
feted chant Mate gration olecences Weed and! iene: numerous 
step size doubling and halving iterations. To prevent 
excessive Lp a eine the tolerances in some cases were 
* Siganaea shedeDy nerodueins eet éorer al “Most. trajectories - 
were relatively unchanged; however, the most sensitive 
output parameters dealt with plenum pressures and motions 
in the vertical plane. Final runs with turns and waves 
used a fixed step size of 0.005 seconds which was determined 
to be the average step size used by the Runge-Kutta-Mersen 
method. At the time of this writing, it seems that more 
accurate results could have been obtained with a 0.001 
second step size while sacrificing four times as much 
computer time. For the purposes of this study, it is felt 


that the .005 second step size was sufficient. 


B. RECOMMENDATIONS FOR FUTURE STUDY 

The first and foremost objective is to determine if the 
loads and motions simulation is valid., When this study was 
initially undertaken it was with the idea of eventually 


comparing simulation data with actual test data. At the 
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time of this writing test data had not been made available. 
Once these data are compared more confidence will be gained 
in the loads and motions simulation thus enabling modifica- 
tions in the simulation for larger craft. Additionally, 
steps should be taken to modify the simulation for the XR-3 
to further develop understanding and confidence in the pro- 
gram. Once these steps have been taken the next obvious 
een is to ibaell hult4thoudandl ton gross ‘weight ‘CAB/SES.” e 
Once it has been proved that the basic program is 
adequate and that it can be modified, a very interesting 
field oF study) would be am optimal cont ron system. Items 
to be condiered for these controls aside fron standard 
rudder and propulsion apparatus would be fan pressure 
control for the seals as well as the plenum, fins, and any 
other type of control surfaces that would enhance the 
stability control that would be required for high speed 


operations. 
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APPENDIX A 


Figure 1. Forty Knot Turns: 
Roll Angle Versus 
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Figure 2. 


Forty Knot Turns: Pitch Angle Versus Time 
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Forty Knot Turns: 
Yaw Angle Versus Time 
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Eighty-five Knot Turns: 
Roll Angle Versus Time 
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Figure 5. 


Eighty-five Knot Turns: 


Pitch Angle Versus Time 
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Figure 6. Eighty-five Knot Turns 
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Figure 7. One Hundred Knot Turns: 
Roll Angle Versus Time 
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Figure 8. One Hundred Knot .Turns: 
Pitch Angle Versus Time 
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Figure 9. One Hundred Knot Turns 
Yaw Angle Versus Time 
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Figure 10. Forty Knot Turns: 
Turn Rate Versus Time 
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Figure 11. Forty Knot Turns: 


Summation of Axial Forces Versus Time 
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Figure 12. Forty Knot Turns: 
Summation of Lateral Forces Versus Time 
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Figure 13. Forty Knot Turns: 
Summation of Vertical Forces Versus Time 
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Figure 14, Forty Knot Turns: : 
Summation of Roll Producing Moments 
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Figure 15. 


Forty Knot Turns: 
Summation of Pitch Producing Moments 
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SUMMATION OF YAW PRODUCING MOMENTS (x1.0E 04 ft-lbf) 


Figure 16. 


Forty Knot Turns: 
Summation of Yaw Producing Moments 
Versus Time 
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Figure 17. 
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Forty Knot Straight Run with Waves: 
Wave Amplitude Versus Time 
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Figure 18. 
dX = 216 ft 
a= 0.5 ft 


Forty Knot Straight Run with Waves: 
Roll Angle Versus Time 


93$ ES 


81 


HA 


seu 


ee Beer 
¢ : 


1 


PITCH ANGLE (x1.0E-02 degrees) 


35 


4 
a 


N10 


Ne 


Figure 19. Forty Knot Straight Run with Waves: 
Pitch Angle Versus Time 
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Forty Knot Straight Run with Waves: 
Summation of Roll Producing Moment 
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Figure 21. 


Forty Knot Straight Run with Waves: 
Summation of Pitch Producing Moments 
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WAVE FORCE X-DIRECTION VERSUS TIME (x1.0E 03 1bf) 
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Figure 22. Forty Knot Turns: Wave Force X-direction 


Versus Time 
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Figure 23. Forty Knot Turn with Waves: Wave Amplitude Versus 
Time 
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Forty Knot Turn Comparison with 
and without Waves: 
Roll Versus Time 
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Figure 25. Forty Knot Turn Comparison with 


and without Waves: 
Pitch Versus Time 
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